GLOBAL SOLUTION TO THE INCOMPRESSIBLE 
FLOW OF LIQUID CRYSTALS 



XIAOLI LI AND DEHUA WANG 

Abstract. The initial boundary value problem for the three-dimensional incompressible 
flow of liquid crystals is considered in a bounded smooth domain. The existence and 
uniqueness is established for both the local strong solution with large initial data and 
the global strong solution with small data. It is also proved that when the strong solution 
exists, a weak solution must be equal to the unique strong solution with the same data. 



1. Introduction 

Liquid crystals are a state of matter that have properties between those of a conventional 
liquid and those of a solid crystal that are optically anisotropic, even when they are at rest. 
In this work, we are interested in a Navier-Stokes type model for incompressible fluids that 
takes into account the crystallinity of the fluid molecules in the three-dimensional case, 
that is, a nematic liquid crystal model, which can be governed by the following nonlinear 
hydrodynamical system (see [T3"l HH] and references therein) : 

y + u^Vu- fj.Au + VP = -XV • (Vd Vd) , (1.1a) 

Q I 

— + u- Vd = 7 (Ad + |Vd| 2 d) , (1.1b) 
at 

V-u = 0, (1.1c) 

where u G R 3 denotes the velocity, d G S 2 (the unit sphere in R 3 ) the unit- vector field 
that represents the macroscopic/continuum molecular orientations, P £ R is the pressure 
(including both the hydrostatic part and the induced elastic part from the orientation 
field) arising from the incompressibility V • u = 0; and they all depend on the spatial 
variable x = (xi, X2, X3) G R 3 and the time variable t > 0. The term AV • (Vd Vd) in 
the stress tensor represents the anisotropic feature of the system. The positive constants 
fj,, A, 7 stand for viscosity, the competition between kinetic energy and potential energy, and 
microscopic elastic relaxation time or the Debroah number for the molecular orientation 
field, respectively. We set these three constants to be one since their exact values do not 
play any role in our analysis. The symbol Vd0 Vd denotes a matrix whose (i,j)-th entry 
is d Xi d ■ d Xj d for 1 < i, j < 3, and it is easy to see that Vd Vd = (Vd) T Vd, where 
(Vd) T denotes the transpose of the 3x3 matrix Vd. 
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System (jl.ip is a simplified version, but still retains most of the interesting mathematical 
properties (without destroying the basic nonlinear structure) of the original Ericksen- 
Leslie model ([3 El QUI CEH EE [14]) for the hydrodynamics of nematic liquid crystals; see 
|15[ [2"0| [22] for more discussions on the relations of the two models. Both the Ericksen- 
Leslie system and the simplified one describe the macroscopic continuum time evolution 
of liquid crystal materials under the influence of both the velocity u and the orientation 
d which can be derived from the averaging/coarse graining of the directions of rod- like 
liquid crystal molecules. In particular, there is a force term in the u-system depending on 
d; the left hand side of the d-system stands for the kinematic transport by the flow field 
while the right hand side represents the internal relaxation due to the elastic energy. In 
many situations, the flow velocity field does disturb the alignment of the molecule, and in 
turn, a change in the alignment will induce velocity. 

We consider the initial-boundary value problem of system (jl.ip in a bounded smooth 
domain Qct 3 with the initial condition: 

(u, d) | t=0 = (u (x), d (x)), xe(], (1.2) 

and the boundary condition: 

(u,^d) |an= (0,0), t>0, (1.3) 

where v is the outer unit-normal vector field on <9S7, uo : £1 — > M 3 , and do : Q — > § 2 are given 
with compatability; for the velocity u the non-slip boundary condition, i.e., homogeneous 
Dirichlet type, is considered, and for the orientation vector d the homogeneous Neumann 
boundary condition is posed here. 

Roughly speaking, system (II. lj) is a coupling between the incompressible Navier-Stokes 
equations and the transported flow of harmonic maps. There have been many studies on 
system ([11]), see [TUl CLU CE3 13 QH Q21 CEJ EH E21 [23] and the references therein. Recently, 
in Lin-Lin- Wang [17], they established both interior and boundary regularity theorem for 
such a system in dimension two under smallness conditions. And, they also established 
the existence of global weak solutions that are smooth away from at most finitely many 
singular times in any bounded smooth domain of M 2 . In Lin-Liu [15] . they addressed both 
the regularity and existence of global weak solutions to the n-dimensional (n = 2, 3) Leslie 
system of variable length, i.e., the Dirichlet energy 

- / |Vd| 2 (fx for d : S n_1 

is replaced by the Ginzburg-Landau energy 

J Q| V d| 2 + (1 ~[2 |2)2 ) dx for d:fi^R". 

More precisely, they proved the global existence of weak solutions with large initial data 
under the assumptions that Uo G L 2 ($7),do G if 1 (fi) with do|an G H?.(d£l) in dimension 
two and three. The existence and uniqueness of global classical solution was also obtained 
if uo G ff 1 (n),do G H 2 (Q) in dimension two or dimension three when the fluid viscos- 
ity \i is large enough. The similar results were obtained also in [22] for a different but 
similar model. When weak solutions are discussed, the partial regularity theorem of the 
weak solution was investigated in [16] (and also [H]), similar to the classical theorem by 
Caffarelli-Kohn-Nirenburgh |3j for the Navier-Stokes equations. 
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In this paper, we are interested in the existence and uniqueness of global strong solution 
(u,d,P) of CULJ in W 2 ' q (nf x W 3 >i(Q) 3 x W^ q {n) with q > 3. By a Strong Solution, 
we mean a triplet (u, d, P) satisfying (jl.lh almost everywhere with the initial-boundary 
conditions f)1.2j) - f)1.3[) . Our strategy to consider (jl.ip is to linearize it as 

^ - Au + VP = -v • Vv - V • ((Vf) T Vf) , (1.4a) 
ot v 

3d 

— - Ad = -v Vf + |Vf| 2 f, (1.4b) 

V-u = 0, (1.4c) 

for some given functions v G M 3 and f 6 M 3 . One of the motivations of making such an 
linearization is that we can use the maximal regularity of Stokes equations (cf. Theroem 
13. 2p and the parabolic equation (cf . Theroem 13. If) . We first use an iteration method to 
establish the local existence and uniqueness of strong solution with general large initial 
data. Then we prove the global existence by establishing some global estimates under 
the condition that the initial data are small in some norm. As system (jl.lj) contains 
the Navier-Stokes equations as a subsystem, one cannot expect generally better results 
than those for the Navier-Stokes equations. The uniqueness of global weak solution is 
still an open problem. We shall prove that when the strong solution exists, all the global 
weak solutions must be equal to the unique strong solution, which is called the weak- 
strong uniqueness. Similar results were obtained by Danchin [3J for the density-dependent 
incompressible Navier-Stokes equations. We shall establish our results in the spirit of [3J, 
while developing new estimates for the orientation field d. Due to the specific structure 
of the equations for u, especially the strongly nonlinear term (Vd) T Ad in the u-system, 
it will be necessary to obtain more regularity for d. 

The rest of the paper is organized as follows. In Section 2, we state our main results 
on local and global existence of strong solution, as well as the weak-strong uniqueness. In 
Section 3, we recall the maximal regularity for Stokes equations and the parabolic equation, 
and also some L°° estimates. In Section 4, we give the proof of the local existence. In 
Section 5, we prove the global existence. Finally in Section 6, we show the weak-strong 
uniqueness. 



2. Main Results 

In this section, we state our main results. If k > is an integer and p > 1, we denote 
by W k,p the set of functions in L p (fi) whose derivatives of up to order k belong to L p (f2). 
For T > and a function space X, denote by L p (0,T; X) the set of Bochner measurable 
X-valued time dependent functions / such that t — > \\f\\x belongs to L p (0,T). Let us 
define the functional spaces in which the existence of solutions is going to be obtained: 

Definition 2.1. For T > and 1 < p,q < oo, we denote by MJ^ q the set of triplets 
(u, d, P) such that 

u G C([0, T];D 1 ~/' P ) n L p (0, T; W 2 «{£1) n W^ q (n)), ^ G 17(0, T; L«(fi)), V • u = 0, 

3(1-1) r),] 

deC([o,r];fC P p, )n^(o,r ; ^(n)), — g if(o,T;L«(n)), 
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and 

P€LP(0,T;W 1,ff (n)), / P dx = 0. 

in 

The corresponding norm is denoted by || • ||m|,> 9 - 

We remark that the condition P <ix = in the definition f|2. 1 1) holds automatically 
if we replace P by 

. . l-ip 

in (jl.ip . Also, in the above definition, the space D A p stands for some fractional domain 

of the Stokes operator in L q (cf. Section 2.3 in [1]). Roughly, the vector-fields of D A p 
are vectors which have 2 — | derivatives in L 9 , are divergence- free, and vanish on <9f2. The 

3(1--) 

Besov space (for definition, see [2]) B QtP p can be regarded as the interpolation space 
between L q and W 3 ' q , that is, 



3(1— i*) 

p 9 ; p p = (l 9 , w r3 ' ,? ) 1 _ i 



p"' 



3(l-i) „ 1 

Moreover, we note that Pg iP p VF 1 ' 9 if p > |. By the imbedding ' 9 L°° as q > 3, 

3(l- i ) 

one has B Q)P p «— )■ L°°, which will be used repeatly in this paper. 

The local existence will be shown by using an iterative method, and if the initial data 
are sufficiently small in some suitable function spaces, the solution is indeed global in time. 
More precisely, our existence result reads: 

Theorem 2.1. Let Q be a bounded smooth domain in M 3 . Assume | < p, q < oo with 



1 - f) G (0, 1) and u G D^**", d 6 bJ£ p) . Then, 

(1) There exists Tq > 0, sitc/i i/tot, system (jl.ip mi/i £/ie initial-boundary conditions 
(|1.2p - (|1.3p /ias a unique strong solution (u, d,P) G Mj.^ 9 infix (0,To); 

(2) Moreover, there exists 5o > 0, such that, if the initial data satisfy 

||Uo|| !-!,„< 5 , ||d || 3(1 _1) < 5 , 

n p fs v 

U A q C 9.P 

i/ien (|1.1|) - (|1.3|) /ias a unique strong solution (u, d,P) G M^' 9 in Q x (0,T) /or aZZ 

r > o. 



According to [17] . a Weak Solution to (jl.ip with the initial-boundary conditions (jl.2p - 
(|1.3p means a triplet (u, d, IT) satisfying system (jl.ip in x (0, T) for < T < +00 in 
the sense of distributions, i.e, for any smooth function ip(t) with i/j(T) = and </>(x) G 
(P^(ft)) 3 with V • cf) = 0, we have 

- / (u,V>V) c &+ / (u- Vu,ip(j))dt + ii I (Vu,ipV(j))dt 
Jo Jo Jo 

= ^(0)(u o ,<£) + A / (Vd Vd, V>V0) (it, 
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and 

r-T f -T f -T 



(d, ip'cj)) dt+ f (u • Vd, dt + 7 [ (Vd, ^V0) dt 
Jo Jo 



= ^(O)(d o ,0)+7 / |Vd| 2 (d,^)dt, 
J o 

where (•,•) denotes the inner product in L 2 (il) 3 . Moreover, (u, d) satisfies ()1 .3|) in the 
sense of trace. 

Next, we will give a uniqueness result. For < T < +oo, suppose (u, d,II) with 
u G L 2 '°°(ft x [0, T\) n W2' (O r ), d G L°°([0, T], i? 1 (O)) n L 2 ([0, T],H 2 (Q)), 

and 

Vn G L5(0,T;LK(n)) 
is a global weak solution to 01)- (fL3l) . Since Vd G L 2 (0,T; H 1 (U)) and |d| = 1, then 

Ad-d + |Vd| 2 = 0. 

Hence |Vd| G L 4 (0 x [0,T]). We have the following energy inequality (cf. |17] . Section 5 
for the two-dimensional case): 



1 [ (|a(t)| 2 + |Vd(t)| 2 )dx+ f [ (|Vu| 2 + |Ad + |Vd| 2 d| 2 )dxds 

2 Jn Jo Jn 

< \ I (|u | 2 + |Vd | 2 )dx, 
1 Jn 



(2.1) 



for all t G (0, oo). We remark that the assumption on pressure function holds since II can 
be determined as in the Navier-Stokes equations (see [9]). 

As for the standard Navier-Stokes equations, the question of uniqueness in the above 
class remains open. However, for the same initial-boundary conditions, a relation between 
the weak solution and the strong solution can be formulated as: 

1 — -,p 3(l-i) 

Theorem 2.2. Assume that uq G D a p and do G B q ^ p p . Then any weak solution to 
(|l.ip - (|1.3p in the above class is unique and indeed is equal to its unique strong solution. 

Usually, we call this kind of uniqueness as Weak-Strong Uniqueness. For the similar 
results on the compressible Navier-Stokes equations, we refer the readers to [61 118|. 



3. Maximal Regularity 

In this section, we recall the maximal regularities for the parabolic operator and the 
Stokes operator, as well as some L°° estimates. 
For T > 0, 1 < p, q < oo, denote 

W(0,T) := (W l > p (0,T;L q {n))f D (LP(0,T;W 3 > q {n))f. 

Throughout this paper, C stands for a generic positive constant. 

We first recall the maximal regularity for the parabolic operator (cf. Theorem 4.10.7 
and Remark 4.10.9 in pQ): 
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Theorem 3.1. Given 1 < p,q < oo, uj q £ B^p p) and f € (L p (0, T; L 9 (IR 3 ))) 3 , the 
Cauchy problem 



9w Aw = /, t > 0, 



) dt 

|w(0) = UJQ 
has a unique solution uj G W(0, T), and 

\\u\\yV(0,T) < C ( \\f\\LP(0,T;Li(R3)) + ll w o|| 3(1-1) J , 

where C is independent of ujq, f and T . Moreover, there exists a positive constant cq 
independent of f and T such that 

IM| W ( 0j t) > c o sup ||w(t)|| 3 (1 _I } . 

te(o,T) p 

Now we recall the maximal regularity for the Stokes equations (cf. Theorem 3.2 in [3]): 
Theorem 3.2. Let Q be a bounded domain with C 2+e boundary in R 3 and 1 < p,q < oo. 
Assume that uq 6 D a J and f G L P (M + ; L q (£l)). Then the system 

^-Au + VP = /, ^Pdx = 0, 
V-u = 0, u| a n = 0, 
u| t=0 = u 

/ias a unique solution (u,P) satisfying the following inequality for all T > 0: 



U A q 



T ' A du\ 



p X p 

dt 

Li 



<C \ \\u \\ Dl _ hp + \\f(t)\\l q dt 



(3.1) 



with C = C(q,p, c(O)) where o~(£l) stands for the open set 
with <5(f2) denoting the diameter of Q. 

Remark 3.1. We notice that (|3.ip does not include the estimate for ||u|| iP (o,T;L<?)- Indeed, 
thanks to u\qq = 0, Poincare's inequality, and the fact J n Vu dx = 0, we have 

1 1 u-l | v^ 2 ,^ < C||Au||x, 9 , 

and then (|3.ip can be rewritten as 

i 

V 



D, 



u(T)|| x i + / ||(VP,u,Au,9iu)||^dt 



<c(l|uo|l Dl -i, P +(^ T ||/(t)||^^ 
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We have the L°° estimate in the spatial variable as follows (cf. Lemma 4.1 

3g 

II jrllS 

LP(0,T;IV 2 .5)' 



Lemma 3.1. iei 1 < p, q, r, s < oo satisfy < | — ^ < 1, £foen 

l|V/|| LP( o,T;L~)<CT5-*r''^ 



L°°(0,T;D A /' ) 

/or some constant C depending only on Q,p,q and 

1-6 1 3 



p 2 2q 

Similarly, we have, 
Lemma 3.2. Let 1 < p, q < oo satisfy < | — || < 1, i/ien 



I|V/||lp(o,T;L°°) < CT 



2 2<j 



L°°(0,T;S,, p p ) 

/or some constant C depending only on £l,p,q and 

1-6 1 3 



LP(0,T;W 2 >4)> 



p 2 2q 

Proof. First, we notice that (cf. Theorem 6.4.5 in [2]) 

]_2_3 x _3 1 — 6 1 3 
(5oc,£ % Boo,^)e,l = SSo.i With = - — — 

p 2 2q 

Also the imbedding 5^ 1 L°° holds due to Theorem 6.2.4 in Hence, 

iv/ioo < ciiv/Hflo t < q|v/|| V 2 liv/ii 1 -^ 

-"00,00 -C>oo,oo 

We remark that (cf. Theorem 6.2.4 and Theorem 6.5.1 in [2]) 

2(1 — — ) 2— — — - 1 — — 

R P <^_v R P 9 W 1 ' 1 ? <^-v R 1 <^_v R 8 



Thus, according to (|3.2p and by applying Holder's inequality, we deduce that 

l|V/|| iP(0 ,T;Loo) <C[f ||V/f e i _3 llV/H^-f 3 dt ) ' 



1 

T \ P 



• I / ILfll^JI/ftl^dt 



< c 



B p q I 



^CT^II/II 1 



LP(0,T;iy 2 ><2)- 
L°°(0,T;B q , p p ) 



The proof is complete. 
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Remark 3.2. In this paper, we will use the following weaker result: 

l|V/|| LP(0 ,T;L-) < CT^WfW 1 - 6 ll/lli^T;^). 

L°°(0,T;B q ) p p >) 

Lemma 3.3. Let 1 < p, q < oo satisfy < § — | < 1, then 



\Af\\LP(0,T;L°°) < CT E 



i i 



3 9 



L°°(0,T;B 9 , P p ) 

/or some constant C depending only on £l,p,q and 

1-8 1 1 



3(i_I) \\J \\LP{0,T;W 3 -i)i 



p 3 q 

Proof. First, we notice that (cf. Theorem 6.4.5 in [2]) 

2_3_3 2—- 1 (9 1 1 

{B^l ',-800,^)0,1 = £^1 with = - - -. 

P ^ q 

Also the imbedding x W 1,0 ° is true due to Theorem 6.2.4 in [2]. Hence, 

II A/IUoc < C||V/|| w i„» < C||V/|| B i < C||V/f 2 _, IIV/HV 3 3 • (3-3) 

D CO t OO -DoO.OO 

We remark that (cf. Theorem 6.2.4 and Theorem 6.5.1 in [2]) 

3 3 r, 3 



3(l-i) 3---- 9 2-'- 

R P' ^ R P 9 W 2 ' 1 ? ^ R 2 <^_v R 9 
r -Doo,oo j r -Dq oo r -Doo,oo- 



Thus, according to (|3.3[) and by applying Holder's inequality, we deduce that 

IIA/H LP(0 ,T;Loo) < C ( f T \\Vf\f 2 3 \\Vff-l~P 3 dt) " 
\Jo B « B p~« / 

^cfrii/cji/ii^^v 

^CT^ill/ll 1 - 6 ll/llLpww"..)- 
The proof is complete. □ 



4. Local Existence 

In this section, we prove the local existence and uniqueness of strong solution in Theo- 
rem 12.11 The proof will be divided into several steps, including constructing the approx- 
imate solutions by iteration, obtaining the uniform estimate, showing the convergence, 
consistency and uniqueness. 
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4.1. Construction of approximate solutions. We initialize the construction of ap- 
proximate solutions by setting 

d° := d , u° := u . 

For given (u n ,d n ,P n ), the Stokes equations (|1.4ap and the parabolic equation (|l,4b|) 
enable us to define (u n+1 , d n+1 , P n+1 ) as the global solution of 



Ssgtl _ Au n+1 + VP n+1 = -u n • Vu n - V • ((Vd n ) T Vd n ), 
dd^± _ Ad n+1 = -u n • Vd n + |Vd n | 2 d n , 



(4.1) 



V-u" +1 = 0, J n P n+1 dx = 
with the initial-boundary conditions: 



n+l jn+l^ 



(u n+i ,d 



|t=0 



(u ,d ), (u n+1 ,cU" +1 )| 9n = (0,0). 



According to Theorems 13. Iti3.2l an argument by induction yields a sequence 

{(u n ,d n ,P n )} neN cM™ 

for all positive T. 



4.2. Uniform estimate for some small fixed time T. We aim at finding a positive 
time T independent of n for which {(u n , d n , P n )} ng pj is uniformly bounded in the space 
M^' q . Indeed, applying Theorem 13.21 to 



9u" 



dt 



Au n+1 + V pn + 1 = _ u n . y u n _ y . ((Vd n ) T Vd n ), 



V-u n+1 = 0, J n P n+1 dx = 0, 

„n+l| „n+l| n 

u T \t=o — u , U \gn — U 



and Theorem 13. II to 



Qd n+1 



Ad 



n+l 



— JU --u n -Vd n + |Vd n |M n , 

d"+ 1 | i=0 = d , cU n+1 U = 0, 



we obtain 



|u" +1 (T)|| 



T 



VP n+1 ,u n+1 ,Au n+1 



<9u 



n+l 



3* 



dt 



Li 



< C I ||u |U_i iP + (jf ||u n • Vu" + V • ((Vd n ) T Vd n ) \\ p Lq dt} 



(4.2) 



and 



d' t+1 (T)|| 3(1 _ 4) + || d ^|| w(0iT) 



n+l I 



B ' p 



< (' I UdolLsd-i) + 11 -u"-Vd n + |Vd n | 2 d n || iP(0)T;L , ) . . 



B 1 p 



(4.3) 
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Now define 



U n (t) : = \\u n \\ j i + U \\LP(0,t;W 2 <o) + \\~aT jy(0,t;£«) 

+ ||d n || 3(^1) + ||d"||w(0,t)j 

L°°(0,t;B q ) p p ') 



and 

U° := ||u || + lldoll 3(i-i) 

Hence, from (|4.2p and (|4.3p . we get, using Lemmas 13.1113,31 



C/ n (£) < C(J7 + ||u n || L oo( 0jt;L9 )||Vu n || L p( 0it;L oo) 

+ 2||Vd n || L oo( 0it . L9 ) ||Ad ||zp(Q,t;Z°°) 

+ ll u?1 ||L° o (0,i;L'J)l|Vd rl || L p( 0jt . L oo) 

+ l|d' 1 ||L° o (0,i;L° o )l|Vd n || iO o( 0it . L 9)||Vd n || L p( 0it . L oo)^ 

< C (V° + 2*3"^ (?7 n (i)) 2 + 2*5~« (C/ n (t)) 2 + ^-25 (t/ n (t)) 3 ) . 
Moreover, if we assume that C/ n (t) < 4Cf7° on [0,T] with 



\ — 



1 < T < 



3 \ 9-3 



(4.5) 



i 64C 2 C/° + 64C 3 (J7 )^ 

then a direct computation yields 

U n+1 {t) < 4CU° on[0,T]. 

From (|4.2p - ()4.4p . we conclude that the sequence {(u n , d n , P n )}^ =1 is uniformly bounded 
in My ? . More precisely, we have 

Lemma 4.1. For all t G [0, T] with T satisfying (j4.5|) . 

U n (t) < 4CU°. 
4.3. Convergence of the approximation sequence. 

Lemma 4.2. {(u 71 , d n , P™)}^ =1 is a Cauchy sequence in and thus converges. 
Proof. Let 

u n := u n+1 - u n , d n := d n+1 - d n , P n := P n+1 - P n . 
Then, the triplet (u n ,d n ,P n ) satisfies 

' ^ - Au n + VP n = -u^ 1 • Vu n - u"- 1 • Vu"- 1 - V • ((Vd ft_1 ) T Vd n ) 

-V^Vd^^Vd"- 1 ), 

— - Ad n = -u 11 " 1 • V 

+ ((Vd n + Vd™- 1 ) : Vd™- 1 ) d 71 " 1 , 
V • u n = 0, J n P n dx = 



2£ - Ad n = -u"" 1 • Vd n - u"- 1 • Vd™- 1 + \Vd n \ 2 d n - 1 (4.6) 
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with the initial-boundary conditions: 

(u n ,d n )\ t=0 = (0,0), (u",^d")|an = (0,0) 

Define 

,du n 

. + Hu'll/xo.tiW 2 ^ + I 



U n (t) :_||u n ||^^ ^ + ll un |lif(0,t;H/ 2 .9) + \\-Q^\\LP(0,t;Li) 



+ \\^ pn \\LP(0,t;L") + lid"!! 3(1-1) + ||d n ||w(0,t)- 

By using Lemmas I3.HI3.3| we obtain the following estimates: 

1 1 — 1 \7„ra i ,,n— 1 V7i»n— 111 

||U • VU +U • VU \\LP(0,t;L<i) 

< ||u" 1 ||L°°(0,t;L9)||Vu n ||ip( 0) t ; Loo) + ||u™ 1 || L oo( 0ji;i q) || Vu™ 1 ||z,P(o,t;L°°) (4.7) 

< 4CU° (^-^llu^lUco^w) + || Vu"- 1 IlLP^ioo)) , 



II v • ((vd^ 1 ) 1 vd") + v • ((vd"- 1 ) 1 vd^)\\ LP{0)t . Lq) 

< ||Vd n || Iy cx! ( - 0)t . L ^||Ad ri_ 11^(0^00) + ||Vd n_1 || LO o( 0j4;L9 )|| Ad n || iP ( 0it;i oo) 

+ ||Vd n_1 || L oo( 0it;L q)|| Ad n_1 || LP ( 0it;L oo) + || Vd n ~ \\l°°(p,p,l«) II Ad n- \\nt>(p t p,L°°) 

< 8CU° (HAd^llip^L-) +ti-«||Vd n - 1 |Ucc (0>tiW) ) , 

llu^-vd^ + u^-vd"- 1 !)^^) 

"—Ill 1 1 V7«JTl II i II, ,1 — 111 IIV7^3"' — 1| 



(4.8) 



^ ll 11 ™ llL°°(0,t;L9)ll^d n || L p( 0i t. L cx)) + ||u n ||L°°(0,t;L'3)l|Vd n \\LP(0,t;L°°) (4.9) 
k°°(0,t;L8) + ||Vd n 1 \\LP{0,t;L°°)) ; 



< 4CC/° t5-35||u n - 1 |Uc O(0>t;W ) + || Vd"" 



and 



U| Vd n| 2a „-i + (( V d n + vd"- 1 ) : vd"- 1 ) d"- 1 !!^^.^) 

< {4CU ) 2 (ts-^Hd^lUoo^rco) + 2||Vd n - 1 || L p (0it;Loo) 
Applying Theorems ELKS to ([46]) . with the help of (j47F]) - ([4T0"j) . we have 

C7"(t) < 4C7C/ (2^-^||u"- 1 || Loo(0>t;M) + || Vu"- 1 !!^^;^) + 2||Ad n - 1 || L p (0it;iO o ) 



(4.10) 



+ 2i3-||Vd"- 1 || L o O(0ii;L<?) + ||Vd"- 1 



L°°(Q,t;Li) "Til llLP(0,t;L°°) 
n-li 



+ 4C7^ u (^-^||d"- i || L o O(0it;Loo) + 2||Vd n - i || L p (0it;Loo) 



(4.11) 



Combining (|4.1ip and Lemmas I3.1M3.31 one has 



U n (t) < 16CU ((1 + 3C77°)P~^ + ts-lj U n - l (t). 

Thus, if we choose To satisfying (|4.5p such that the condition 

16C77°(2 + 3CU°)T ^~« < -, or 16C77°(2 + 3C7[/°)T 5-4 < _ 

is fulfilled, it is clear that {(u n , d n , P n )}^ =1 is a Cauchy sequence in and thus con- 
verges in ■ □ 
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4.4. The limit is a solution. Since {(u n , d n , P n )}^ =1 is a Cauchy sequence in M^ q , 
then it converges. Let (u, d, P) e M^ q be the limit of the sequence {(u n , d n , P n )}^_ 1 in 
Mj? 9 . We claim all those nonlinear terms in f)4. 1[) converge to their corresponding terms 

in ([Lip in (L p (0, Tq; L^fJ))) 3 . Indeed, due to the convergence of u n to u in M™, we have 

||u n - Vu n -u- Vu|| LP(0iTo;L9) 

= || (u n - u) • Vu" + u • V(u n - u)|| l p (0jTo;L9) 

< ||u n — u|| i oo( 0jTo; z / <j)||Vu n || i p( 0jTo;i oo) + ||u|| L oo( 0iTo;Z/9 )||Vu n — Vu|| iP ( 0jTo;i0 o) 

< 4CC/°T 5~4 || u n - u\\ m p,« + CT^~h ||u||£c~ ( o,To;W)l|u n - u|| M ™ 
-> 0. 

Hence, 

u n -Vu n ^u-Vu in (L p {0,T ;L q {n))) 3 . 

Similarly, we have 

V • ((Vd") T Vd") -> V • ((Vd) T Vd) in (L p (0, T ; L 9 (^))) 3 ; 

u n Vd"^uVd in (L p (0,T ;L q (n))) 3 . 

Since 

|||Vd"| 2 d"-|Vd| 2 d|| LP(0jTo . £8) 

< |||Vdf d™ - |Vd"| 2 d|| L p (0 ,T 0; W) + |||Vd"| 2 d - |Vd| 2 d|| L p (0 ,T ;M) 

< ||d n — d|| i0 o( 0jTo;i oo- ) ||Vd ?1 || i o C (o i T'o;L9)||Vd r! '|| L p(o i T ;L° ) 

+ l|d||L° o (o,T ;L°°)||Vd re + Vd|| L o ( 0jTo;L q)||Vd ri - Vd|| LP ( 0iTo;L oo) 

< {(4cu ) 2 + ||d|| LOO(0iTo;LOO) (4Cc/° + ||Vd|| LOC(0iTo;W) )) rt^lK - d\\ M ™ 

-+0, 

then 

|Vd n | 2 d n -> |Vd| 2 d in (L p (0,T ;L q (n))f. 

Thus, taking the limit as n — > oo in f|4.1|) . we conclude that (jl.ip holds in L p (0, T ;L q (Q)), 
and hence almost everywhere in $7 x (0,Tq). 

Multiply the d-system, i.e., (jl.lbp by d, we obtain 

^ + ^u.V(|d| 2 ) = Ad.d + |Vd| 2 |d| 2 . 

Since 

A(|d| 2 ) = 2|Vd| 2 + 2d-(Ad), 

then it follows that 

l^f + l U • V(|d|2) = ^ A(|d|2) " |Vd|2 + |Vd|2|d|2 - 
Therefore, it is easy to deduce that 

d (|d ^~ ^ - A(|d| 2 - 1) + u • V(|d| 2 - 1) - 2|Vd| 2 (|d| 2 - 1) = 0. (4.12) 
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Multiplying (|4.12p by (|d| 2 — 1) and then integrating over fi, using f 1 1 . 1 c [) and (jl.3p . we 
get the following inequality: 



4- I (|d| 2 -l) 2 dx<4 / |Vd| 2 (|d| 2 - l) 2 dx 

< 4||Vd|||oo f (|d| 2 - l) 2 dx. 
Jn 



(4.13) 



Remark that interpolating between L°°(0, T ; W^ q ) and L p (0,T ;W 3 > q ) shows that for 
some positive f3 > \, d belongs to T 2 (0, T ; H 2+/3 ) and that llVd^ 6 L^O,^). Notice 
that 

/ (|d| 2 - l) 2 dx = 0, at time t = 0. 
Jn 

Thus, using (|4.13p together with GronwaU's inequality, it yields |d| = 1 in O x (0,T ). 

4.5. Uniqueness. Let (ui,di,Pi) and (u 2 ,d 2 ,P 2 ) be two solutions to (11.11) with the 
initial-boundary conditions (|1.2|) - (|1.3|) . Denote 

u = ui - u 2 , d = di-d 2 , P = Pi - P 2 . 

Note that the triplet (u, d, P) satisfies the following system: 

( g - Au + VP = -u • Vui - u 2 • Vu - V ■ ((Vdi) T Vd) - V • ((Vd) T Vd 2 ) , 
§ - Ad = -ui • Vd - u • Vd 2 + |Vdi| 2 d + ((Vdi + Vd 2 ) : Vd) d 2 , 
Vu = 0, / Q Pdx = 

with the initial-boundary conditions: 

(u,d)| t=0 = (0,0), (uAd)|an = (0,0). 

Define 

, . .. _ ,. .. _ .. .. du 

X(t) := U j_l + U \\LP(Q,t;W 2 'i) + ^7 LP(0,t;L«) 

L°°(0,t;D A(j p ) K " ' Ot 



+ \WP\\Lv(0,t;Li) + l|d|| 3(l-i) + ll d ll>V(0,i)- 

Thus, repeating the arguments in (|4.7p - (|4.10p . we have 

X(t) < WCU° ((1 + 3CU°)t*~% + i* - *) 

< ix W . 

Hence, X(t) = for all t £ [0, To], which guarantees the uniqueness on the time interval 
[0,T ]. 
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5. Global Existence 

In this section, we prove that, if the initial data is sufficiently small, the local solution 
established in the previous section is indeed global in time. To this end, we first denote 
by T* the maximal time of existence for (u, d, P). Define the function H(t) as 

. . .. .. .. .. ii^ u n 

H(t) := U !_i + U \\ LP r ,t ;W 2 >1) + "ST \\LP(0,t ;L") 

L°°(0,t;D Aq ? P ) 8t 

+ \\^P\\LP(0,t;L9) + 3(i-I) + ||d||w(0,t)- 

L°°(Q,t;B q ; p p ') 

To extend the local solution, we need to control the maximal time T* only in terms of 
the initial data. For this purpose, it is obvious to observe that H(t) is an increasing and 
continuous function in [0, T*), and for all t S [0,T*), we have, using Theorems [XTE21 

(5.1) 

(5.2) 

(5.3) 
(5.4) 

(5.5) 



H(t) <C(U° + ||u • Vu|| iP(0ii;L<?) + || V • ((Vd) 1 Vd) || L P (0 , t ;L,) 
+ ||u ■ Vd|| iP ( 0it;i g) + || I Vd| 2 d|| iP ( 0it;i9 )^ . 
On the other hand, Lemmas 13.1113.31 imply that 

||U • Vu|| L p( 0j4;L g) < ||u|| L oo( 0it;i q)||Vu|| L p( 0it;L oo) 

|| V- ((Vd) T Vd) \\Lp(o,t;Li) < C||Vd|| £ , O o( 0) t ;£ , 9 )||Ad|| £ ,p( 0)f . £ ,oo) 

||U • Vd|| L p( 0j4;L g) < ||u|| L oo( 0it;i q)||Vd|| L p( 0it;L oo) 

and for the fact that |d| = 1, we have 

|||Vd| 2 d|| iP (o i4; L9) < ||Vd||xo (Q )t . i ?)||Vd|| i p(o,t;L 00 ) 

Substituting (f5T2]) - (f575]) into (^T|) . we get 

H(t) <c(l7° + (3^"4 + t^«)H 2 (t)^ . (5.6) 
Assume that T is the smallest number such that 

H(T) = ACU°. 

This is possible because H{t) is an increasing and continuous function in time. Then, 

H{t) < H(T) = 4CU°, for all t E [0, T), 
and from (|5.6p . we deduce that 

16C 2 U° ( 3T^"^ + rW) > 3. 



This implies that the maximal time of existence T* will go to infinity when the initial 
data approaches zero. More precisely, we can show that, if the initial data is sufficiently 
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small, the solution exists globally in time. To this end, we need some other estimates for 
the terms on the right side of (|5.ip . Indeed, by W 1,q L°° as q > 3, we have 

||U • Vu|| L p( 0j4;L g) < ||u|| L oo( 0it;i ,)||Vu|| L p( 0jt;Lo o) 

< CH 2 {t). 

Similarly, 

||V- ((Vd) T Vd)|| LP(0it;i9) , ||u- Vd\\ L v { Q,tm, IHVd| 2 d|| iP(0ii;Lq) < CH 2 {t). 
Thus, (|5.ip turns out to be 

H(t)<C(U° + H 2 (t)). (5.7) 
Now we take U° sufficiently small such that 

Then, we compute directly from (|5.7|) and the continuity of H(t) that 



„, , 1 - VI - 4C 2 U° 1 

4,(4)5 2C S 2C 

for all t 6 [0,T*), which implies ||(u, d, P)\\M v ' q bounded. Hence, according to the local 
existence in the previous section, we can extend the solution on [0, T*) to some larger 
interval [0, T* + rj) with r\ > 0. This is impossible since T* is already the maximal time of 
existence. Hence, when the initial data are sufficiently small, the strong solution is indeed 
global in time. 

The proof of Theorem 12.11 is complete. 



6. Weak-Strong Uniqueness 

The purpose of this section is to show Weak-Strong Uniqueness in Theorem 12.21 To 
this end, we first formally deduce and obtain an energy estimate for the strong solution 
to jHHES). 

Lemma 6.1. Letp,q satisfy the same conditions as Theorem \2.1\ and (u, d, P) £ MjS g be 
the unique solution to (jl.ip on 0, x [0, T] . Then for any < t < T , we have 

(|u(t)| 2 + |Vd(t)| 2 )dx + 2 f [ (|Vu| 2 + |Ad+ |Vd| 2 d| 2 )dxcis 
i Jo Jn 



[ (|u | 2 + |Vd | 2 )dx. 
Jn 



Proof. Note that 

u G C([0, T); D A ^ P ) n L p (0, T; W 2 ' q n Wq 



d € C(%T]-B^ ^)nLP(0,T;W^), 
D A ~ V,V ^ B^p~» n X q (see Proposition 2.5 in @]), 

where 

X q = {z£ L q (nf | V • z = in ft and z • n = on <9ft}, 
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for some positive /3 > |, it follows from the standard interpolation inequalities that 
u G C([0, T];HP) n L 2 (0, T; H 1+ ?), u G L 4 (fl x [0, T]), 

d G C{[0,T};H 1+ ' 3 )r\L 2 (0,T;H 2+ls ), Vd G L 4 (ft x [0,T]). 

This enables us to justify the following computations. 
Multiplying (jl . la|) by u, integrating over J7, we get 

~ f \u\ 2 dx + f \Vu\ 2 dx = - f u • ((Vd) T Ad) dx. (6.1) 
2 at Jq Jq Jq 

Here we have used the facts 

V • (Vd Vd) = V P^Y^) + (Vd) T Ad, 
and V • u = in S7, u = on d£l, as well as 

u • Vu • u dx = [ VP u dx= [ V ^ • u dx = 0. 



2 

Multiplying (|l.lbj) by —(Ad + |Vd| 2 d) and integrating over $7, we obtain 

- ( ^ • Ad dx - [ (u • Vd) • Ad dx = - [ | Ad + |Vd| 2 d| 2 cix. 
Jn dt Jq Jq 

Here we have used the fact that |d| = 1 to get 

dd , „ V7a\ IV7^|2j _ 1 ( IV7J|2^l d | 2 , „ V7lJ|2|V7J|2 



(6.2) 



+ u • Vd • iVdrd = - I Vdr^r-!- + u • Vldrl Vdr = 0. 
dt J 2 V dt J 

Since d u d = on d£l, integrating by parts, we have 

[^■Addx = -l± [ \Vd\ 2 dx. 
Jn dt 2dtJ n ]l 

Hence we obtain 

~ [ \Vd\ 2 dx+ ( |Ad+ |Vd| 2 d| 2 dx= ( (u • Vd) • Ad dx. 

2 dt Jq Jq Jq 

By adding (|6.ip and (j6.2|) . we eventually get the identity: 

14- I d u l 2 + |Vd| 2 )dx + / (|Vu| 2 + | Ad + |Vd| 2 d| 2 )(ix = 0, (6.3) 

2 dt Jq Jq 
for all t G (0,T]. 

Integrating f)6.3f) over the time interval [0, t], we obtain the energy equality of this 

lemma. □ 

We remark that ()6. 3f) is usually called the basic energy law governing the system (jl.lj) - 
(|1.3j) . It reflects the energy dissipation property of the flow of liquid crystals. 
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Now we proceed to prove Weak-Strong Uniqueness. In view of the regularity of u, we 
deduce from the weak formulation of (jl.la[) that 



u • u dx + / / Vu : Vu dxds 
< > Jo Jn 

t r ft 



[ |u | 2 dx- f [ (Vd) T Ad • u dxds + [ [ u • + u • Vu) dxds. 
Jn Jo Jn Jo Jn V vs ) 



On the other hand, since u satisfies (jl . la|) . i.e., 



^ = Au - u • Vu - VP - V ) - I Vd) Ad. 



then we have, from (|6T 



ii'|u'Vu+ (Vd) T Ad - u • Vu) dxds. 
Similarly, in view of the regularity of d, we have 



Vd : Vd cbc - / |Vd | 2 dx 
<> Jn 

t r 

d • Ad s + u ■ Vd • Ad - Ad • Ad - |Vd| 2 d • Ad) duds. 

o Jn 

Taking advantage of (jl.lbp . we obtain, from (|6.6 



o Jn 



(6.4) 



u • u dx — / |uo| 2 dx 
n Jn 

2 Vu -.Vu dxds -II (Vd) T Ad • u dxds (6.5) 
Jo Jn Jo Jn 

io Jn 



(6.6) 



Vd : Vd dx - / |Vd | dx 
n Jn 

t r , 

2Ad- Ad + u- Vd- Ad + u- Vd- Ad ( 6 - 7 ) 
- |Vd| 2 d • Ad - |Vd| 2 d • l\d \ dxds. 



From (|2.ip . (j6.5[) . (|6.7p and the fact that (u, d,P) (resp. (u, d,IT)) is a strong solution 
(resp. weak solution) to (jl.ip with the initial-boundary conditions ()1.2j) - (|1.3|) .we have the 
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following energy estimate of (u — u, d — d): 

\ [ (\u-u\ 2 + \Vd- Vd| 2 ) dx 
* Jn 

<- f /(|Vu- Vu| 2 + | Ad- Ad| 2 ) dxds 
Jo Jn 



rt 

io Jn 




(u - u) ■ Vu ■ (u - u) + (Vd - Vd) ■ Ad ■ (u - u) 

- u • (Vd - Vd) • (Ad - Ad) 

+ (|Vd| 2 d - |Vd| 2 d) • (Ad - Ad)) dxds 

(|Vu - Vu| 2 + | Ad - Ad| 2 ) dxds + /, 



o Jn 



where 



t 

' (u - u) • Vu • (u - u) + (Vd - Vd) -Ad-(u-u) 




io Jn 



- u • (Vd - Vd) • (Ad - Ad) 
+ (|Vd| 2 d - |Vd| 2 d) • (Ad - Ad)) dxds. 
Next, we will estimate I term by term. By the zero boundary condition, we have 

— / (u — u) ■ Vu ■ (u — u) dx = / (u — u) ■ Vu • u dx 

Jn Jn 

= / (u - u) • (Vu - Vu) • u dx 

Jn 

< ||u||x,oo ||Vu — Vu|| L 2 1| u — u|| L 2 



1 _ llu|| 2 

< 2 ll Vu - Vu lli^ + ^^ll u - u lli^ 



- / (Vd - Vd) • Ad • (u - u) dx < || Ad|| L -> ||u - u|| L 2 ||Vd - Vd|| L 2 
Jn 



S 

Jn 



<^fc(||u-a|| 2 2 + ||vd-vd||i 2 ), 

u • (Vd - Vd) • (Ad - Ad) dx < ||u|| L oo || Ad - Ad|| L 2 ||Vd - Vd|| L 2 

<|l|Ad- Ad|| 2 2 + M^||Vd - Vd|| 2 2, 



- f (|Vd| 2 d - |Vd| 2 d) • (Ad - Ad) dx 
Jn 

< ||Vd|||oo||d-d|| L2 ||Ad- Ad|| L a 

+ || Vd + Vd|| L oo || Vd - Vd|| L 2 1| Ad - Ad|| L 2 

< l\\Ad - Ad|| 2 2 + ||Vd|||oo||d - d|| 2 2 + ||Vd + Vd|||oc ||Vd - Vd|| 2 2 . 
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Then, we eventually get from (|6.9p - (|6,12p that 

/ < / ' ( i||Vu - Vu||f 2 + || Ad - Ad||| 2 + IHI '°° + ii Ad il^ || U _ fi" 2 

H 

~2 



L 2 



+(||vd + vdHi. + l|u|li °° + ' |Ad||L °° ) || va - vd||| 2 ( 6 - 13 ) 



+||Vd||i 0o ||d-d||i a jcfa. 

Now, we wish to estimate ||d — d||^2 . We write 

d t (d - d) + u • V(d - d) + (u - u) ■ Vd = Ad - Ad + |Vd| 2 d - |Vd| 2 d. (6.14) 
Multiply (16.141) by d — d and integrate over 0, x (0, t), we have 

- / |d-d| 2 dx 

= -/ / (u-u) • Vd- (d-d) dxds- / / |Vd - Vd| 2 dxds 
Jo Jn Jo Jn 

+ [[ |Vd| 2 |d - d| 2 dxds + [ [ (Vd + Vd) : (Vd - Vd) d • (d - d) d*ds. 
Jo Jn Jo Jn 

Using Sobolev's inequality ||u — u||^6 < C|| Vu — Vu||^2 and for some e > small enough, 
it is easy to get 

1 



2 



|d-d| 2 dx 

n 

t 



< [ fc e ||Vd|| 2 3 + ||Vd|||oo + [ |d-d| 2 dxds + e [ [ |Vu-Vu| 2 dxds 
Jo V 2y J n J J n 



(6.15) 



2 

< 



Now we have from (f6T8|) (f6TT3j) and (pJ5|) that 

1 f (|u - a| 2 + |d - d)| 2 + |vd - vd)| 2 ) 

2 Jn 

t ll u ll|°° + l|Ad|| Loo f 2 

/ / |u — u| axds 

Jo 2 Jn 

+ J i\ + II Vd Hi» + C » Vd »!» + ll vd lli") J l d " d l 2 d * d > (6.16) 

+ r»} + Mj, + ||Ad||^ + 3||vd + vdirj,\ i |yd _ ya|2 dxds 



< C [ ip(a) [ (|u - ti| 2 + |d - d| 2 + |Vd - Vd| 2 ) dxds, 
Jo Jn 

where 

<p{s) = 1 + ||Vd||ioo + ||Vd||loo + ||Vd||i» + ||u||i» + ||Ad|| L ~. 
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Notice that HVdH^oo, HuH 2 ^,, ||Ad||i<x> £ L 1 (0,t). Moreover, by applying the quasi- 
linear equations of parabolic type estimates (cf. [12J Chapter VI, Section 2) to (jl.lbp . we 
see d(-,i), d(-,i) G C 1,a with respect to the space variables, for some a > 0, and its C 1,a 
norm is independent of t. Then we have f(s) S L 1 (0,t). Applying Gronwall's inequality 
to (|6.16p . we obtain 

f (|u - u| 2 + |d - d| 2 + |Vd - Vd| 2 ) (fx = 
Jn 

for all t. Thus, u = u, d = d a.e. and P = II up to a constant in SI x (0, T). 
The proof of Theorem 12.21 is now complete. 
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